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We present a chronological review of the progress in calculating weak matrix elements using staggered fermions. 
We review the perturbative calculation of one-loop matching formula including both current-current diagrams and 
penguin diagrams using improved staggered fermions. We also present preliminary results of weak matrix elements 
relevant to CP violation calculated using the improved (HYP (II)) staggered fermions. Since the complete set of 
matching coefficients at the one-loop level became available recently, we have constructed lattice operators with 
all the g 2 corrections included. The main results include both AI = 3/2 and AI — 1/2 contributions. 



1. INTRODUCTION 

Staggered fermions preserve enough chiral sym- 
metry to prevent those operators of our physical 
interest from mixing with operators of wrong chi- 
rality and to protect quark mass from additive 
renormalization. This chiral symmetry is essen- 
tial to calculate e'/e. Staggered fermions have the 
advantages over DWF of requiring less comput- 
ing time and dynamical simulations are already 
possible with relatively light quark masses. 

By construction, staggered fermions possess 
four degenerate flavors (also called "tastes"). 
They allow large taste changing quark-gluon in- 
teractions, which make the perturbative correc- 
tion large even at the one loop level. In addi- 
tion, four-fermion operator mixing matrix is huge 
(65536 x 65536), which makes it impractical to de- 
termine the matching coefficients in a completely 
non-perturbative way. Another disadvantage is 
that the unimproved staggered fermion action 
and operators receive large perturbative correc- 
tions |1I2| and have large scaling violations of or- 
der a 2 [SJ- Both of these disadvantages can be 
alleviated by improving staggered fermions using 
fat links. It has been shown that taste symme- 
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try breaking in the pion spectrum is significantly 
reduced with fat links |4I5| . 

The first goal of this staggered e'/e project is 
to check the results (e'/e < 0) obtained using 
quenched domain wall fermions (DWF) by CP- 
PACS and RBC collaboration 7 . The sec- 
ond goal is to extend the calculation to dynami- 
cal simulation. The main goal is to find a window 
for new physics or to confirm the standard model, 
when our numerical results are compared with the 
experimental results. 

This project has progressed through the follow- 
ing small steps. 

2. PERTURBATION 

Since the penguin diagrams had already been 
known for the gauge invariant operators from 
|5] , we calculated the remaining diagrams of the 
current-current type £Q. This provided a com- 
plete set of matching formula for e'/e, combined 
with the existing results for penguin diagrams [S] . 
We also observed that the perturbative correc- 
tion for some of the operators such as O7, 0% 
and Oq are so large that the uncertainty from the 
truncated contribution is of order 100% for using 
unimproved staggered fermions. 

3. NUMERICAL STUDY 

Using unimproved staggered fermions, we per- 
formed a numerical study on e'/e. Using the 



matching formula given in |1I8| . we constructed 
fully one loop matched gauge invariant operators. 
As expected, we found large perturbative correc- 
tions for and B§, so that we can not quote 
quantitative results for this although the statisti- 
cal uncertainties are under control. We also ob- 
served that different quenching transcriptions of 
the continuum operators on the lattice (proposed 
by Golterman and Pallante [5] lead to noticeably 
different values for Bq ^U] . This indicates a large 
quenching uncertainty, which deserves further in- 
vestigation. 

4. IMPROVEMENT SCHEME 

Hence, it was essential to reduce the large per- 
turbative correction by improving the staggered 
fermions. Hence, the main goal was to find an 
improvement scheme which can reduce the per- 
turbative correction down to 10% or less. In or- 
der to find the best scheme, we calculated, ex- 
plicitly, one loop matching coefficients for var- 
ious improved staggered actions and operators 
P]: 1) Fat7, 2) Fat7+Lepage, 3) HYP and 4) 
AsqTad-like (Fat7+Lepage+Naik) actions. We 
observe that all the above improvement schemes 
significantly reduce the size of matching coeffi- 
cients. After the second level of mean-field im- 
provement, the HYP and Fat7 links lead to the 
smallest one-loop corrections Since the HYP 
fat link reduces the non-perturbative flavor sym- 
metry breaking more efficiently in the pion spec- 
trum than the Fat7 link, we adopted the HYP 
scheme in our numerical study. 



5. HYP/Fat7 

We studied further on the HYP link [T2]. The 
HYP links possess some universal properties, 
which are summarized in the following 5 theo- 
rems. 

Theorem 1 (SU(3) Projection) 

Any fat link can be expanded in powers of gauge 
fields (A^). 

B^ =0{A n ) 



1. The linear term, B^ is invariant under 
SU(3) projection. 

(2) 

2. The quadratic term, is antisymmetric 
in gauge fields. 

Theorem 2 (Triviality of renormalization) 

1. At one loop level, only the B^' term 
contributes to the renormalization of the 
gauge-invariant staggered fermion opera- 
tors. 

2. At one loop level, the contribution from 
B^ for any n > 2 vanishes. 

3. At one loop level, the renormalization of the 
gauge-invariant staggered operators can be 
done by simply replacing the propagator of 
the A^ field by that of the B^ ] held. 

This theorem is true, regardless of details of the 
smearing transformation. 

Theorem 3 (Multiple SU(3) projections) 

1. The linear gauge field term Bjp in the per- 
turbative expansion is universal. 

2. In general, the quadratic terms may be dif- 
ferent from one another. But all of them 
are antisymmetric in gauge Relds. 

3. This theorem is true, regardless of the de- 
tails of smearing. 

Theorem 4 (Uniqueness) 

If we impose the perturbative improvement condi- 
tion of removing the flavor changing interactions 
on the HYP action, the HYP link satisfies the 
following: 

1. The linear term Bp in perturbative expan- 
sion is identical to that of the SU(3) pro- 
jected Fat 7 links. 

(2) 

2. The quadratic term B^ is antisymmetric 
in gauge helds. 

Theorem 5 (Equivalence at one loop) 

If we impose the perturbative improvement con- 
dition to remove the flavor changing interactions, 
at one loop level, 



Operators [P x P][P x P] n Operators [S x P][S x P] n 

NAIVE 2 x (111.3 -2CW) NAIVE 2 x (95.6 - 6C N ) 

HYP/Fat7 2 x (6.925 - 2C H ) HYP/Fat7 2 x (19.1 - 6C H ) 



Table 1 

One-loop correction to (03)7/. 

1. the renormalization of the gauge invariant 
staggered operators is identical between the 
HYP staggered action and those improved 
staggered actions made of the SU(3) pro- 
jected Fat7 links (often called W'j, 

2. the contribution to the one-loop renormal- 
ization can be obtained by simply replacing 
the propagator of by that of . 

The first two theorems were used in [2] and 
[T5| . although they did not present their deriva- 
tion. For derivations of all five theorems and fur- 
ther details, refer to 12 . As a result of these 
theorems, we can prove that for each Feynman 
diagram, 

II C fat || < || c thm 11 . (1) 

Here, Cf a t (Cthin) represents perturbative cor- 
rections to gauge-invariant staggered operators 
constructed using SU(3) projected fat links (thin 
links). This inequality is not valid for those fat 
links without SU(3) projection. Hence, this lead 
to a conclusion that we may view the SU(3) pro- 
jection of fat links as a tool of tadpole improve- 
ment for the staggered fermion doublers ■ We 
also present alternative choices of constructing fat 
links to improve staggered fermions in ^21- The 
above five theorems make the perturbative calcu- 
lation simpler for the HYP scheme, because one 
can perform the calculation merely by replacing 
the thin link propagator with that of the HYP 
links. This simplicity is extensively used in calcu- 
lating the renormalization constants of the four- 
fermion operators in the next stage. 

6. PERTURBATION FOR HYP/Fat7 (1) 

We calculated the current-current diagrams to 
obtain perturbative matching coefficients for the 
staggered four-fermion operators constructed us- 
ing the HYP/Fat7 links ^31- m particular, we are 



Table 2 

One-loop correction to (03)7/. 

interested in the (O3)// operator (we use the same 
notation as in U), because this receives large per- 
turbative corrections (~ 1) in the case of unim- 
proved staggered fermions. 

(Oa)// - 2([P xP][PxP]-[Sx P][S x P]) n 

In Tables Q] and |3 we present values of one- 
loop corrections to (O3)// constructed using both 
the HYP/Fat7 link and unimproved thin link 
(denoted as NAIVE). Here, C N = 13.159 and 
Ch = 1.4051, which correspond to tadpole im- 
provement contributions. The results shows that, 
by choosing the HYP/Fat7 scheme, the pertur- 
bative corrections are reduced to « 10% level. 
In the case of the HYP/Fat7 scheme, note that 
the size of one-loop correction is already under 
control even without tadpole improvement. For 
further details, refer to |14) . 

7. PERTURBATION FOR HYP/Fat7 (2) 

In order to obtain a complete set of matching 
formula to convert the lattice results into con- 
tinuum values, we need to calculate the penguin 
diagrams for the gauge- invariant staggered opera- 
tors constructed using various fat links jTH] . One 
of the main results regarding the diagonal mixing 
from penguin diagrams can be summarized into 
the following theorem: 
Theorem 1 (Equivalence) 
At the one loop level, the diagonal mixing co- 
efficients of penguin diagrams are identical be- 
tween (a) the unimproved (naive) staggered op- 
erators constructed using the thin links and (b) 
the improved staggered operators constructed us- 
ing the fat links such as HYP (I), HYP (II), Fat7, 
Fat7+Lepage, and Fat7. 2 

2 Note that AsqTad is NOT included on the list. In this 
case, by construction the operators are made of the fat 
links which are not the same as those used in the action 
due to the Naik term. In addition, the choice of the fat 
links are open and not unique. 
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Figure 1. Bk{h = 1/a) 



The details on the proof of this theorem will be 

given m puma. 

By construction, gluons carrying a momen- 
tum close to k ~ 7r / 'a are physical in staggered 
fermions and lead to taste changing interactions, 
which is a pure lattice artifact. In the case of 
unimproved staggered fermions, it is allowed to 
mix with wrong taste 1) and the mixing coef- 
ficient is substantial. In contrast, in the case of 
improved staggered fermions using fat links of our 
interest such as Fat7, Fat7 and HYP (II), the off- 
diagonal mixing with wrong taste vanishes and 
is absent. In the case of the improvement using 
HYP (I) and Fat7 + Lepage, the off-diagonal mix- 
ing with wrong taste is significantly suppressed. 
The details of this off-diagonal mixing will be 
given in |16| . 



8. NUMERICAL STUDY FOR HYP/Fat7 

Recently we have performed a numerical study 
using Columbia QCDSP supercomputer with the 
HYP (II) staggered fermion. Since we know the 
complete set of matching formula for HYP/Fat7 
staggered operators from |14I15| . we constructed 
fully one-loop matched gauge-invariant operators 
to study e'/e. Here, we present preliminary esti- 
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Figure 2. sf J - 3/2 (^ = 1/a) 



mates oiB K , Bf ,2 \ B (3/2) and b£ /2) calculated 
using the HYP (II) staggered fermions at (3 = 6.0 
on a 16 3 x 64 lattice with 218 configurations. 

8.1. B K 

Fig. n snows Bk as a function of M K , where 
the mesons are made of degenerate quarks. We 
fit Bk to the form suggested by quenched chiral 
perturbation theory J7]: B K = c + Ci(M K ) 2 + 
C2(Mk) 2 log(Mft-) 2 . The cross symbol in Fig. 
1 corresponds to the value interpolated to the 
physical kaon mass. Our preliminary result is 
Bk = 0.581(18), which is consistent with the 
continuum extrapolated value calculated using 
unimproved staggered fermions .'! . However, our 
value is substantially different from the value 
(Bk — 0.6790(16)) calculated using unimproved 
staggered fermions at finite lattice spacing (1/a = 
2.01 GeV). Therefore, this leads to the prelim- 
inary conclusion that, using the HYP staggered 
fermions, the scaling behavior is so improved that 
our results at j3 — 6.0 (1/a = 1.95 GeV) is al- 
ready in agreement with the continuum extrapo- 
lated values of the unimproved staggered results. 
Of course, this claim needs further investigation 
at weaker couplings. In the chiral limit, we obtain 
cq = 0.13(15), which is also consistent with those 
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Figure 3. B^- 3/2 (^ = I /a) 



Figure 4. B^ 1 ' 3 ' 2 ^ = I /a) 



results obtained using the NLO, large N c calcu- 
lation [18119] . The value for Cijc§ is consistent 
with the predictions of quenched chiral perturba- 
tion theory |17| within large errors. 



8.2. b\ 3/2) and B {3/2) 

A major contribution to the AI = 3/2 am- 
plitudes comes from Bj 3 ^ 2 ^ and B^ 2 \ We fit 



-By 3 ^ 2 ' and Bg >/Z> to the form of our trial function: 



,(3/2) 



Kj -Co + Cl {M 2 K ) + c 2 (M 2 K ) \og{M 2 K ). The 



results of^ 3/2) 



and Bg are presented in Fig.0 
and Fig. 01 respectively. Our preliminary values 
at the physical kaon mass are By 1 ' = 0.919(13) 
1.047(15). We also fit B (3 /2) and 



and B 



(3/2) _ 



(3 /2) 

B s to the form suggested by the quenched 
chiral perturbation h{M%) = c + a(M%). 
These fitting results of Bf^ 2 ^ and B^ 2 ^ are 
shown in Fig. and Fig. [51 respectively. Our pre- 
liminary values at the physical kaon mass from 
these fits are b\* /2) = 0.917(13) and B (3/2) = 
1.044(15). Hence, we observe not only that the 
interpolated values at the physical kaon mass are 
insensitive to the fitting functions, but also that 
the extrapolated values in the chiral limit are ex- 
tremely sensitive to the fitting functions. We can 
think of two possibilities for the chiral extrapola- 



s 



tion. One possibility is that the data point at the 
lightest quark mass might be shifted due to a fi- 
nite volume effect, which certainly need further 
investigation in near future. Another possibil- 
ity is that the truncated higher order corrections 
from the quenched chiral perturbation theory are 
not negligible so that we might have to include 
higher order terms to do better chiral extrapo- 
lation. For these reasons, we do not quote our 
values for B^ 2 ^ and B^ 2 ^ in the chiral limit 
but we only quote the values interpolated to the 
physical kaon mass as above. 

(3 /2) 

Note that we calculated B\ at the scale 
ji = X/a using the HYP staggered fermions, which 
would not have been meaningful for unimproved 
staggered fermions due to large perturbative cor- 
rections. Compared with previous calculation 
done using Landau-gauge operators [23j, the sys- 
tematics of the HYP staggered operators are sig- 
nificantly reduced and the results are as much 
more reliable. 

8.3. B^ /2) 

A major contribution to AI — 1/2 ampli- 



tudes comes from B, 



(1/2) 



There are two inde- 

(1/2) 



pendent methods to calculate B v & 'in (partially) 
quenched QCD: the standard (STD) method 
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and the Golterman-Pallante (GP) method UJ. 
First, we fit Bg 1/2 ^ to the form suggested by 
the quenched chiral perturbation theory 
h(M 2 K ) = c + Cl {M 2 K ) + c 2 {M 2 K ) log(M|). The 
results are presented in Fig. [5] Our preliminary 
values at the physical kaon mass are 

B^ 1/2) ' STD (fi= l/a) = 0.714(91) 
B [ 6 1/2) ' GP (n = l/a) = 0.974(69). 

where l/a = 1.95 GeV, set by the p meson mass. 
We also fit f?g 1//2 ' ) to another form suggested by 
the quenched chiral perturbation UJ: ^(Afjj-) = 
c + Ci (Afg-). The results are presented in Fig. [7| 
Our preliminary values at the physical kaon mass 
from this fit are 

B£ /2) ' STD (n=l/a) = 0.701(93) 
B { e /2hGP {^= l/a) = 0.973(70). 

As in the case of B^ 2 ^ and B^ 2 \ we observe 
that the extrapolated values in the chiral limit are 
extremely sensitive to the fitting function choice. 
At present, we can not exclude the possibility that 
the data point at the lightest quark mass could be 
shifted due to a finite volume effect (in the e re- 
gion of the chiral perturbation theory) . It is also 
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Figure 6. B^ 1 ' 1 ' 2 ^ = l/a) 



possible that truncated higher order corrections 
from the quenched chiral perturbation theory are 
so significant that we might have to include these 
neglected terms for chiral extrapolation. There- 

1 /2 

fore, we cannot quote our values of B 6 ' in the 
chiral limit. 

Note that unlike the unimproved staggered 
fermion calculations where the perturbative cor- 
rections to the STD calculation are rj 50%, the 
perturbative corrections in this calculation are 
modest for both the STD and GP methods. We 
also observe that the gap between the STD and 
GP methods is reduced at the physical kaon mass 
using the HYP staggered fermions compared that 
of the unimproved staggered fermions. 

8.4. Preliminary e'/e 

We use the formula provided in [21] to convert 
(m s +m d ) n , B^ 2 \^ = m c ) and B^ 2 \„ = 
m c ) at m c = 1.3 GeV into e'/e. As mentioned 
before, we have two different methods to cal- 
culate B£ /2) in quenched QCD: the STD and 
GP methods. When we use the STD method 
for ^ 1/2 V = m c ), we obtain e'/e(STD) = 
0.00046(23). For the GP method for B^ 1/2) (^ = 
m c ), we obtain e'/e(GP) = 0.00115(17). These 
values are very preliminary and we have not in- 
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eluded an analysis of the systematic errors. In ad- 
dition, we did not use lattice values for any B 



(1/2) 



except for B\ 
them. 



(1/2) 



since we have not yet extracted 



9. FUTURE PERSPECTIVES 

We plan to calculate all the B^ 2 ^ and incorpo- 
rate all of them into the calculation of e'/e. We 
also plan to obtain the optimal matching scale, 
q* [22] • We plan to extend our calculation to dy- 
namical simulation. Current calculation of e'/e 
uses the leading order chiral perturbation results 
to convert K — > tt and K — ► amplitudes into 
K — » 7T7T amplitudes. In principle, it is possible to 
calculate K — > tttt amplitudes directly on the lat- 
tice. A number of attempts in this direction have 
been tried |5S 24 . There has been a numerical 
study on Fat7 as a new fat link computationally 
cheaper than the HYP link, while preserving all 
the nice features of the HYP link |25| . 
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